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Abstract 

We develop an analytical method for studying the properties of a non-interacting Wormlike Chain 
(WLC) in confined geometries. The mean field-like theory replaces the rigid constraints of confinement 
with average constraints, thus allowing us to develop a tractable method for treating a WLC wrapped 
on the surface of a sphere, and fully encapsulated within it. The efficacy of the theory is established 
by reproducing the exact correlation functions for a WLC confined to the surface of a sphere. In 
addition, the coefficients in the free energy are exactly calculated. We also describe the behavior of 
a surface-confined chain under external tension that is relevant for single molecule experiments on 
histone-DNA complexes. The force-extension curves display spatial oscillations, and the extension of 
the chain, whose maximum value is bounded by the sphere diameter, scales as / _1 at large forces, in 
contrast to the unconfined chain that approaches the contour length as J -1 / 2 . A WLC encapsulated 
in a sphere, that is relevant for the study of the viral encapsulation of DNA, can also be treated 
using the MF approach. The predictions of the theory for various correlation functions are in excellent 
agreement with Langevin simulations. We find that strongly confined chains are highly structured 
by examining the correlations using a local winding axis. The predicted pressure of the system is in 
excellent agreement with simulations but, as is known, is significantly lower than the pressures seen 
for DNA packaged in viral capsids. 

1 Introduction 

The Wormlike Chain (WLC) model [1], that well describes the elasticity of DNA, microtubules, and 
polyelectrolytes, and is suitable for a polymer with two length scales: the contour length L and 
persistence length l p . When confined to the surface or volume of a sphere, a third length scale appears: 
R, the radius of confinement. The emergence of this new length scale drastically alters the behavior of 
the WLC, by restricting the conformational space available to the polymer. A strongly confined WLC 
will adopt a tightly bent configuration, which is energetically unfavorable in bulk conditions. Because 
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a number of biologically relevant systems involve stiff chains in confined geometries or adsorbed onto 
curved surfaces, a general understanding of the WLC model in these geometries is essential. 

The confinement of biopolymers to a curved surface is of interest in a number of systems. In 
cukaryotes, the first level of chromosomal compaction of DNA (with l p w 50nm) is histone wrapping, 
with the DNA wrapped around the cylindrical histone (with radius 4.2nm and height 2.4nm) [2, 3, 
4]. The stability of the tightly bent structure is essential in understanding the development of the 
chromosome. In addition, many authors have studied the behavior of polymers confined to the surface 
of both cylinders [5, 6, 7] and spheres [5, 8, 9, 10] to discern the free energy scaling and equilibrium 
behavior of surface confined chains. In particular, an exact solution for the end-to-end distance R of 
a WLC confined to the surface of a sphere has been determined [8], and confirmed using simulations 
[9]. 

Experiments on the dsDNA-containing bacteriophages ^29 [11, 12] and el5 [13], as well as the T 
[14, 15, 16, 17, 18], P [19], and A [20, 21, 22, 23] phage classes, have determined a number of details of 
the structures, pressures, and ejection timescales of many viruses. Regardless of the shape of the viral 
capsid, it is generally seen that the DNA orders itself in concentric rings [14, 20, 13, 15, 16, 21, 24], with 
the spacing <~ 0.3nm between rings. Single molecule experiments [25] have shown the pressure on the 
capsid walls to be on the order of 60 atm, inducing a significant resistance to the DNA encapsulation. 
These observations have generated a number of theoretical studies [26], primarily interested in the 
packaged structure [27, 28, 29, 30, 31], inter-strand spacing [32, 33, 34, 35], energy or pressure [36, 32, 
27, 30, 37, 38, 39, 40], and the loading or ejection process [33, 41, 27, 28, 42, 30]. While the specific 
geometry of the confining viral capsid varies from phage to phage, the properties of the encapsulated 
DNA can be studied using spherical [29, 42, 28] or cylindrical [32, 34] confinement to a very good 
approximation. The study of confined WLCs in these simple geometries is relevant to our understanding 
of the properties of viruses. Interactions between monomers play a significant role in the energetics of 
structure formation in viral packaging. However, the initial stages of the encapsulation process may 
be understood by examining the non-interacting chain, where self-intersections are relatively rare and 
short range interactions simply renormalize the persistence length [43]. 

In order to study the effects of both surface and volume confinement on the behavior of a WLC, 
we will extend the Mean Field (MF) method [44, 45, 46] introduced by Ha and Thirumalai. The MF 
method has been successful in producing tractable theories involving WLCs in many different potentials 
[45, 47, 48]. In particular, it has been applied to the study of a long, closed WLC on the surface of 
a sphere [10]. The paper is organized as follows. In Sec. 2, we extend the MF theory for a WLC on 
the surface of a sphere of radius R, and show that it reproduces all known averages and scaling laws. 
We also show that the theory accurately reproduces the correct scaling coefficient of the free energy 
of confinement. In Sec. 3, we adapt the MF theory to study the behavior of a surface-confined WLC 
subject to an external mechanical force. The application of the MF theory to a WLC encapsulated in 
a sphere (referred to as volume confinement) is discussed in Sec. 4. We show that the volume confined 
chain is in excellent agreement with simulations using the theory. As suspected in previous studies, 
the analytic calculations explicitly show the pressure due to confinement of a non-interacting WLC 
can not reproduce the large values observed in experiments. We also show that the structural order of 
the confined WLC, that is absent in the bulk, can be understood by using a local winding axis. The 
ordering of the chain is purely a consequence of entropic confinement. 
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2 Confinement to the Surface of a Sphere 



Theoretical Considerations 



We begin by developing the Mean Field (MF) formalism for a WLC with persistence length l p , fixed 
inter-monomer spacing a, and length L = Na, confined to the surface of the sphere of radius R. We 
define r„ = (x n ,y n , z n ) the position of the n th monomer, and the bond spacing u„ = Ar„ = r„ +1 — r„ 
with |u„| = a. (shown in Fig. 1). The distribution of the chain in phase space is 

*s({r n }) oc HS(r 2 n -R 2 )6(Arl-a 2 )e- l ^ 3 ^^- A ^ 2 (1) 
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where the second line follows from the first after a Fourier transform of the delta functions (with 
the Fourier variables {A„} and {k n }). Following Ha and Thirumalai [44, 46], we write the partition 
function as Z — j J7 n d 3 r„^/s({r„}) = J ]J n d\ n dk n exp(— Ts[{\ n , k n }]), which defines the free energy 
functional for surface confinement, Ts- The free energy can be written as Ts = T x + + J 7 ? — 
a Y,n(^ri + k n ), where T x is given by = Jl\ n dx„ exp(-H x [{x n }}), with 

n ^ ' 

We assume Ts is sharply peaked around a particular set of Fourier variables {A„,fc„} = {A*,fc*}, 
so that Z ~ exp(— JFg) (i.e. a saddle point approximation). The optimal values of A and k are 
determined by minimizing Ts, i.e. by solving dTs/d\ n — dTs/dk n — 0. In this approximation, the 
Fourier variables play the role of spring constants restricting the position ({k n }) and bending ({A„}) 
of the chain. 

Since the discrete Hamiltonian is quadratic in the x n % we can write Ts = 3/2 log[Dct(Q)] — 
a^2 n (X n + k n ) + const, where the symmetric, (N + 1) x (N + 1) tridiagonal matrix Q is given in 
Appendix A (Eqs. 40-42). The solution for these coupled equations is intractable for large N, and 
additional approximations are necessary to make further progress. The symmetry of the matrix is 
respected by the substitution A„ — > A and k n — > k, except for exactly three elements near the endpoints 
(see Appendix Afor more details). This is similar to the excess endpoint fluctuation terms found in 
the unconfined theory [44, 46], where A was shown to be constant except at the endpoints. With these 
observations, we take 

fci = fcjv+i = k H - 

a R 

k 2 = k N = (3) 

^1 — Ajv — AH — -, 

a R z 

with k n = k and A„ = A for all other values of n. The specific forms of the endpoint terms in Eq. 3 are 
chosen to ensure convergence of the continuum limit. Substitution of these values into the Hamiltonian 
in Eq. 2 and taking the continuum limit (with a — > 0, — > oo, and Na — > L), we can separate the 
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Hamiltonian into interior and endpoint terms, H x = Hq +H e , with 



Ho = £d S (^x 2 ( S ) + Xx 2 ( S ) + k^) (4) 



H e = 5{ul+ul) +11 (^+ X ^+2 l2 (u X ^-u L X ^y (5) 

where we have defined uo — x(0) and ul = x(L), with x = dx(s)/ds. The free energy functional T x 
in the continuum limit becomes 



T x = - log 



J rf 4 xcxp(-W e ) J V[x(s)]eM-Ho[x(s)]) 



(6) 



with x = (x ,xl,uo,ul), and the total free energy is 

Ts = T x + T y +T z -\L-kL-25-2 ll . (7) 
The path integral in Eq. 6 can be evaluated exactly [49], and we find 



Z (x) = J V[x{s)\ exp(-H [x(s)}) = K exp ^x • Mx 



(8) 

where M is a 4x4 matrix; M and K are evaluated in Appendix B, and given explicitly in Eqs. 47 and 
53 in terms of the two frequencies 



Expressions resulting from the propagator in Eq. 8 can greatly simplified in the limit of large Lw i; 
which we refer to as strong confinement (see below). 
The total free energy functional finally becomes 

T s - -3 log ^ J d 4 x Z(x)J — XL — kL — 25 — 2 7l 

Z(x) = Z (x)exp(-W e ) (10) 

with H e given in Eq. 5 and Zq in Eq. 8. The optimal parameters A, k, 8, 71, and 72 are obtained by 
solving the five coupled Mean Field equations, 

d\ dk dS <9 7 i ^72 
Note that, from Eqs. 7 and 11, the A, k, 5 and -fi derivatives immediately imply, respectively, 

L .,„ 

<uS + ui> = 2, (rg + r|)=2i? 2 . 



This suggests that the MF approximation is equivalent to replacing the local requirements u 2 (s) = 1 
and r 2 (s) = R 2 by the global conditions (u 2 (s)) = 1 and (r 2 (s)) = R 2 . The parameter A plays the role 
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of a spring constant that keeps the bond spacing fixed on average, while k is a spring constant that 
keeps (r 2 ) = R 2 on average. The 72 derivative in Eq. 11 implies (ro • uo) — • Ul) = 0, as is expected 
since the rigid constraints require u(s) to be tangential to the surface of the sphere (i.e. u(s) _L f for 
all s). 

The solutions to the mean field equations (Eq. 11) can be determined exactly for all L, l p , and R, 
giving 



A 



9 

8L 



lp 



k = 



lp 

2R 2 ' 



3 
4' 



7i 



3 
4' 



and 72 = -i|. 



(12) 



We note that A changes sign for R 2 < 8l 2 /9, because the bonds tend to be more compressed with 
decreasing R, which requires a net repulsion between neighboring monomers to satisfy the constraint 
(u 2 ) = 1. With the solutions to the MF equations in Eq. 12, the frequencies in Eq. 9 become 



w» = 



3 

4L 



1 ± A/1 - 



16§ 
9R 2 



(13) 



We note that, in the limit of large R, Lw\ ~ L/l p and Loj 2 ~ Ll p /R 2 . Our demarcation of strong 
confinement, Luui ^> 1, requires long chains (L » l p ) and sufficiently small radii (R <C ^Ll p ). 
Correlation Functions 

The bending correlation function can be computed directly using the solutions in Eq. 12. However, 
in the limit as R — ► 00, we find (u(0) ■ u(L)) — * er'' L l 2l '> as R — » oo, rather than the expected two- 
dimensional correlation function, (u(0) • u(L)) = e r L ' 2l p. This suggests that the theory requires a 
mean field persistence length, l , much like in the unconfined theory [44, 46]. Substitution of l p = 3/ 
into the correlation function results in the expected limit as R — > 00. Ha and Thirumalai, who 
found a similar result for a three dimensional unconfined WLC with l p = ZIq/2, argued that the 
renormalization of l p in the MF theory arises because of the additional forbidden chain conformations 
allowed by replacing the 5 functions in Eq. 1 with Gaussians. Consequently, the mean field persistence 
length is smaller than the true persistence length. In the confined theory, we allow three dimensional 
configurations by replacing the confining 8 functions with Gaussians, which would be forbidden by 
the surface confinement, in addition to relaxing the rigid inter-monomer constraints. For this reason, 
we would expect the confined MF theory to permit additional conformations of the WLC that are 
forbidden by the rigid constraints, relative to the unconfined theory, thus increasing l p /lo- In practice, 
l p is often determined by fitting experimental or simulation data to a suitable polymer model. Hence, 
the renormalization of l p within the mean field theory is not a serious concern. 

The correlation functions computed using the MF theory can be written as 



(r(s)-r(s')) = i? 2 e~ |As|/c 
(u(s)-u(s')) - e H As l/fe 



cosh 



|A£| 

Cs 



l~ sinh f 



lAal 



As I 

< 1 isli ( '— — Qs 



V Cs 



n s 



(14) 



with fls = \/l — 161q/R 2 , and the decay length of the correlations (s = l/4Zo. We note that Eq. 
14 has reproduced the exact calculation of Spakowitz and Wang [8], valid for all values of L, lo, and 
R. The ability to calculate these averages exactly shows the accuracy of the MF method. We can 
also verify directly that (u 2 (s)) = (r 2 (s))/R 2 = 1. Higher order moments are incorrect, though, since 
(u 4 ) = (r 4 )/R 4 = 5/3 ^ 1 as the rigid constraints would require. 



5 



The Free Energy of Confinement 

We can determine the free energy of confinement for the system (which does not require the 
substitution of l p = 3lo) as 

PF ~ Ts = W p + 2% + const - (15) 

This result is identical to the scaling predicted by Odijk for a tightly bent WLC [5]. Additionally 
the coefficient of the scaling law agrees with that predicted by Mondescu and Muthukumar [7] for the 
surface-confined Freely Jointed Chain. We use the Configurational Bias Monte Carlo (CBMC) method 
[50] to determine the scaling coefficient of the free energy for L = 50a, for various values of l p and R. 
The theoretical curves for [3F <~ Ll p /2R 2 +const are accurate to within ~5% (see Fig. 2). 



3 Surface Confined Stiff Chains under Tension 

Theoretical Considerations 

The efficacy of the mean field method is in its ability to study the effect of additional potentials in 
problems involving confined WLC's with relative ease. In this section, we apply an external tension, 
f , to the ends of a surface confined WLC. For the free chain, the Mean Field method has been shown 
to give excellent agreement with experimental results [51]. Such a calculation for the surface-confined 
WLC will also verify that the mean field method satisfies the confinement on average, even under the 
extreme situation of a strong pulling force. 

The distribution in phase space of a confined WLC under tension can be written as ^(f) = 
^f s exp[— (3i ■ (ri — r/v+i)], with (3 = 1/ksT and ^5 given in Eq. 1. Because the external tension does 
not generate an energetic term quadratic in the r„'s, the MF theory in the previous section can be 
used with little change. The discrete free energy functional can be written as ^Fs(f) — F x + T y + T z — 
Pf(%N+i — x\) — a^2 n (X n + k n ), with T x given in Eq. 2 (T y and T z are similarly defined), and where 
we have taken f = fx. None of the terms involving k n or A„ are altered with the application of the 
force, and we can again rewrite the quadratic terms of the Hamiltonian using a symmetric, tridiagonal 
matrix Q (explicitly given in Appendix A, Eqs. 40-42). This again suggests the replacement used in 
Eq. 3, with A„ and k n constant except near the endpoints. In the continuum limit (N — > 00, a — > 0, 
and Na — > L) , we find the free energy 



D AL+fcL+2<5+27i 



(16) 



with Z(x) given in eq 10 and x = (xo, Xl, Uq, ul), a result similar to Eq. 10. The integrals can be 
evaluated with little difficulty, yielding 



Fs(f) = ^s(0) + 



(iW) 2 



M 



- M" 1 

11 V / 12. 



(17) 



where ^"s(0) is the free energy at / = (Eq. 10), and M is given in Eq. 47. 

Under the assumption that Luii 3> 1 (a strongly confined chain, see Eq. 9), it is not difficult to 
show that the solutions to the mean field equations (Eq. 11) become 



A 



8L 



l p 



k = 



lp 



7i = 1 



72 = - 



2R' 



(18) 
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Under the application of a force, only the the endpoints of a strongly confined chain arc affected, 
reflected in the fact that only 71 depends on /. The interior monomer behavior should be relatively 
insensitive to / far from the endpoints, so it is not surprising that A and k are independent of the 
force. 

Force- Extension Curves 

The extension as a function of the external tension can be computed using (R) = —dFs / d((3f). 
For a strongly confined chain, we find 



(xl - x ) 
R 



1 



Rfif 



(R 2 )q 

R? 



1 + 



(R 2 )o 

4R 2 



v/1 + 4(iW)2/9 



(19) 



with (R 2 )o = 2R 2 (1 — (r • i"l)o) the average end-to-end distance with / = (with (r • tl)o given in 
eq. 14. While the force-extension curves for a confined WLC increase monotonically as a function of 
/, the system has rather complicated behavior as a function of R. In Fig. 3(a), we see the extension 
of a stiff chain (l p — 150a, approximately the persistence length of DNA) as a function of R is highly 
oscillatory for small R, due to the non-monotonic behavior of (R 2 )o as a function of R. However, 
oscillations in (R 2 )o are not observed in more flexible chains, as seen in 3(b). This is due to the fact 
that (r -yi)q <~ e~ L / 4 ' , so that the oscillations in (R 2 )o for longer or more flexible chains are damped 
out. The non-monotonic behavior observed in 3(a) is thus due to finite-size effects. 
The asymptotic limits of Eq. 19 are 



(X L - Xq) 



/3/(R 2 )o/3 



2R 



&(1 



(R 2 



-) + o(.r 2 ) 



(20) 



In the low force regime, the system has the expected linear response to the tension, and (ii-io) < 2R 
for all values of the force. Surprisingly, though, the scaling of {xl — xo) — 2R <~ / _1 in the high force 
regime differs from the scaling of the unconfined chain, (x^ — x ) — L ~ / _1//2 . The change in the 
large-force scaling laws is linked to the fact that only the endpoints are affected by the force for a 
surface confined chain. For a free WLC, the extension of the chain comes about by alignment of all 
bonds with the force axis. When confined to the surface of the sphere, the extension occurs primarily 
by to the translation of the endpoints to the poles of the sphere, rather than a global realignment of 
the bond vectors. This is reflected in the fact that 71, which controls the position of the endpoints, is 
the only mean field variable dependent on /. We note as well that the / _1 ^ 2 scaling is seen in the MF 
theory for the free WLC, and comes about due to the fact that A (which determines the behavior of 
all of the bonds) becomes a function of /. 

It is also possible to numerically solve the mean field equations for small L/R and R(3f ^> 1 (strong 
stretching limit), where we find (x L — x ) m 2Rsin(L/2R) for L < wR, the exact end-to-end distance 
of a fully stretched chain confined to the surface of a sphere. The mean field method thus satisfies 
the confining constraints on an average, even under high forces, and again predicts the lower moments 
exactly. 

Finally, we can determine the free energy of a confined WLC under tension, in the limit of strong 
confinement: 



(3F- 



9L LI 



2# + 31 ° S 



4(/3/i?) 2 



-3L/4l t 



f3f(x L - x Q ) + const, 



(21) 



The force-dependent terms in Eq. 21 are not extensive because the tension only strongly effects the 
endpoints of the chain. If we neglect terms on the order of e~ 3L/>4L p, F becomes tension-dominated 
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when / exceeds a critical force R/3f c <~ Ll p /AR 2 . Because the only force-scale in the problem is /3PA ~ 
Llp/R 2 , with ^4 the surface area of the sphere, the scaling of this critical force is expected. We expect 
the leading coefficient to be correct, due to the accuracy of our expression for the free energy of a WLC 
without the external tension (see Fig. 2). 

It is amusing to estimate f c for a strand of DNA wrapped around a histone [2, 3, 4], with l p ss 50nm, 
L w 43nm, and R w 4nm. We find the tension dominates the free energy when / > f c w 34pN, which is 
significantly larger than the force required at each unwrapping event seen in single molecule experiments 
on histones [3]. However, as it has been observed that the tilting of the histone with respect to the 
force axis is of great importance when determining the behavior of the system [2], which the mean field 
theory does not take into account. Our result only provides an upper bound on the unravelling force. 

4 Wormlike chains confined to the interior of a sphere 

Theoretical Considerations: 

The mean field theory for computing the average properties of a surface confined chain can be 
extended to studying the effects of volume confinement. The distribution in phase space of a WLC 
confined to the interior of a sphere is 

*y({r„}) oc Yl B(R 2 - v 2 n ) 6{Ar 2 n - a 2 ) e -*p(Ar„ +1 -Ar„) 2 /a 3 (22) 

n 

where Q(x) is the Heaviside step function, that ensures that each monomer is contained within the 
sphere. The last two terms in Eq. 22 are identical to the ones in ^5 (Eq. 1). The similarities between 
the two distributions suggest that volume confinement can be treated at the mean field level as well. 
Unfortunately, the O function in Eq. 22, that ensures the chain is within the interior of the sphere of 
radius R, can not be dealt with as simply as the S functions found in Eq. 1 at the mean field level. 
It is not difficult to show that, for a single particle confined within a sphere, simply minimizing the 
Fourier Transform of the function does not give the correct value of (r 2 ). However, we may formally 
write 

/ioo 
JJ dk n d\ n exp 
-ioo „ 



where g, an undetermined function, is chosen such that free energy minimization satisfies the rigid, local 
constraints on average (u 2 (s) = a 2 and r 2 (s) < R 2 ). We immediately see that the same substitution 
of interior (i.e. k n — > k and A„ — > A) and cndpoint terms (Eq. 3) will satisfy the symmetry of Q 
(see Appendix A, and Eqs. 40-42), due to the similarities between Eq. 1 and Eq. 23. This allows 
the problem of volume confinement in the continuum limit to be written in terms of the mean field 
variables A, fc, S, 71, and 72, with the free energy expressible as Ty = T x + T y -\-T z — G[X, k, S, 71, 72]. 
T x is defined in Eq. 6, and G constrains the minimization of T (i.e. contains the as yet undetermined 
Lagrange multipliers). 

The treatment of volume confinement at the mean field level is more difficult than the case of 
surface confinement for a number of reasons. In the case of surface confinement, we replaced the strict 



1 , (Ar„+i 
-air, 



Ar„) 2 



Ar 2 



R 2 



g(ak n ) 



(23) 
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constraint of r 2 (s) = R with the global constraint j J„ L ds(r 2 (s)) = R 2 . While the average monomer 
position for a volume confined WLC is not known a priori, we expect that interior monomers, those 
far from the endpoints, will have a uniform behavior. This suggests that we write 



d S (r 2 ( S )) = P R\ (24) 



for some unknown p, which may depend on L, l p , and R. Eq. 24 is equivalent to the requirement at 
the mean field level 

% = L P (25) 

(see Eq. 4). Unlike the surface case, the average position of the endpoints within the sphere need 
not be identical to the average position for interior points of the chain, i.e. (r 2 (s)) ^ const for volume 
confinement. At the mean field level, this can be treated approximately by the restriction 

(r 2 ) = (rl)= Po R 2 , (26) 

with po ^ p an unknown parameter. This is implemented using 

wr 2po (27) 

(see Eq. 5). In addition to the nonuniformity at the endpoints, volume confinement allows for different 
fluctuations in the bending at the endpoints. Because u(s) need not be perpendicular to f (as was the 
case for surface confinement), (u(s) • r(s)) ^ 0, which must be accounted for at the mean field level 
as well. Since u(s) = dr(s) /ds changes sign under the transform s — > L — s, it is simple to show that 
(u • r ) = — (ui • r^). In particular, if the endpoints of a confined chain are found near the wall of the 
sphere, the direction of the bond vectors at the endpoint will be restricted, pointing away from the 
wall of the sphere and giving (yl ■ Ul) > 0. We then restrict 

(ui • t l ) = -<uo • r ) = Pc R. (28) 

The unknown parameter p c represents the average correlation between the position and the bending 
at the endpoints of the chain. Eq. 28 is expressed at the mean field level as 

(see Eq. 5). The lagrangc multipliers for both A and 5 remain unchanged at the mean field level, 
with dT '/d\ — L and dT '/dS — 2 (see Eq. 10). The three Mean Field parameters p, po and p c can not 
be computed within the framework of the MF theory, and must be supplied using some other method. 
We use low friction Langevin dynamics simulations to determine the equilibrium behavior of a WLC 
confined to the interior of a sphere. The details of our simulations are given in Appendix C. 

Up to the three undetermined parameters (p, po, and p c ), we can write (as we did in Eq. 7) the 
mean field free energy for volume confinement as 

T v = T x +T y + T z -XL-25- pkL - 2 Poll + A Pcl2 (30) 
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where T x is identical to the one dimensional free energy functional for the surface case (eq 6). The 
mean field equations for the volume case, given in Eq. 11, are quite similar to the surface equations, 
and can be solved in the limit of strong confinement (i.e. Luji 3> 1, see Eq. 9). We find 

_9_ lp_ = l p = 3(po + pi) 

8l p pR? 2p 2 R 2 4( Po -p2) 

_^ | 3 _ lp ip C /„-. \ 

Ap + 2{p -pl) 72 2pR + 2{p -piy [6 > 

and the frequencies in Eq. 9 become LUi = 3/4Z p (l ± ^/l — l6lp/9pR 2 ). The solutions for the Wj's, 

which define the average behavior over the entire length of the chain, are identical to those found for a 
wormlike chain confined to the surface of a sphere of radius ^fp R (see Eq. 13). However, the endpoint 
terms differ from the surface confined system (Eq. 12), allowing for differing behavior between the 
monomers at the ends and those the interior of the chain. 

In Figure 4, we show the simulated average monomer positions as a function of s for varying R and 
l p . Fig. 4a shows that "interior" monomer behavior (where (r 2 (s)) ~ pR 2 = const) begins to emerge 
in the range 2R > s > L — 2R. Significant deviations from (r 2 (s)) w pR 2 occur near the endpoints 
of the chain over a range of s w 2R, due to the differing fluctuations in the endpoint monomers (see 
Eq. 31). The range of the endpoint effect makes physical sense: if r is near the boundary of the 
sphere, the bending energy near the endpoint will be lower if u is directed towards the center of the 
sphere, as opposed to being directed towards the wall. This suggests that segments of the chain near 
the endpoints will be directed inwards, giving rise to the decrease in (r 2 ) seen in Fig. 4. Endpoint 
effects will dominate the behavior of the chain until the segment comes into contact with the opposite 
side of the sphere on average, a distance of at most 2R. 

In Fig. 4b, we see that increasing the persistence length of the chain while keeping R fixed changes 
the values of p and po (reflected in the overall increase in (r 2 )), but docs not significantly alter the 
qualitative behavior of (r 2 ) as a function of s. For strong confinement, fluctuations in (r 2 ) are small far 
from the endpoints, clearly indicative of an effective surface confinement on a sphere of radius y/pR, 
consistent with the results of the MF theory. 

Determination of p, po, and p c : 

We perform a number of simulations in order to determine the mean field parameters p, po and p c 
for varying l p and R, shown in Fig. 5(a). We find that, for long chains (L = 100a and L = 200a, there 
is virtually no variation in any of the mean field parameters with respect to L. Since p determines 
the effective surface confinement (sec Fig. 4b), increasing L does not change p, as the chain simply 
wraps further around the effective surface at ^fp R. p and p c are likewise independent of the length 
of the chain, due to the fact that the behavior of the endpoints is only weakly dependent on L as long 
as L > 2R. Since the only remaining length scales in the system are l p and R, we expect that all 
of the mean field parameters depend only on the ratio l p /R. This is confirmed in Fig. 5(a), as the 
computed values for p, po and p c each collapse onto a single curve as a function of l p /R. We find for 
long, stiff chains (with L/R 3> 1 and l p /R ^> 1) that p < 0.9, p < 0.95, and p c < 0.25. Most systems 
of biological interest (the viral packing of DNA, for example) are in the strongly confined regime. 

Correlation Functions: 

For interior monomers (where the system is confined approximately to a sphere of radius ^fp R, 
see Fig. 4), we find that the bending correlation function converges on (u(s) • u(s')) — ► e - 3 l A;s l/ 2i P 
in the limit of R — ► oo, rather than the expected unconfined limit of (u(s) ■ u(.s')) = e~l As l // ' p . This 
suggests the Mean Field persistence length l p = i/2l a for large R, identical to the result found in the 
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unconfined MF theory [44, 46]. However, for l p /R ^> 1, the system is effectively confined to a sphere 
of radius ^fp R. As the surface MF persistence length is given by l p = 31$ , we expect the volume 
confined l to be a function of R. The ratio l p /lo = a should be independent of L for long chains, 
since we have seen that the development of near-surface confinement depends only on the ratio l p /R 
(Fig. 5a). We expect 3/2 < a(l p /R) < 3, i.e. a wormlikc chain confined to the interior of a sphere will 
behave somewhere in between a free wormlike chain (3 dimensional), and a surface confined wormlikc 
chain (two dimensional). For interior monomers, we find 



(r(s)-r(s')) « pR 2 e-\ As Vt 
(u(s)-u(s')) ~ e - |As|/Cv 



coshf^vV^sinhf^O 



, AA \ 1 /|Aa| 
cosh I — — S2y — — — sinh — — S2y 
Qv / "v V Cv 



(32) 



with £y = 4a/o/3 and f2y = -^/l — 16a 2 ZQ/9pi? 2 . Near the endpoints, the correlation functions become 
more complicated, due to the dependence of the behavior of the endpoints on p and p c . 

In order to determine l p /lo — a(l p /R), we turn to our simulation results again. The simulated 
bending correlation function is fit using Eq. 32, with a as the only a fitting parameter, with the 
resulting values are shown in Fig. 5b. We find that a does indeed vary with only l p /R for weak 
confinement, and a > 3/2 in this range. For stronger confinement, the bending correlation function 
is only weakly dependent on a, with large fluctuations in the fitting parameter for increasing l p /R. 
However, the saturating value appears to be a(oo) w 5/2. The fact that a does not reach the maximal 
value of lo/l p — 3 is not surprising, as volume confinement still allows fluctuations in (r 2 ) forbidden 
by surface confinement. The decay length in Eq. 32, Qy = 4aZ /3 ~ 3.3Z for strongly confined 
chains, is strictly less than than the decay length for surface confinement, (s = 4Zo, again due to the 
larger number of configurations that are available to volume confined chains. The agreement between 
simulation and theory is excellent not only for the bending correlation function (Fig. 6), but also the 
agreement for the position correlation function is equally as good (data not shown). 

Probes of structures: 

It is of interest to probe the confinement-induced structure in a WLC, that is a thermally fluctuating 
filament in the bulk. Information about the structure of a stiff chain confined to the interior of a sphere 
can be determined using the local winding axis of the chain. The unit local winding axis of bonds i 
and i + 1 (the axis about which and u, + i wind) is given by [9] sn — aj/|aj|, with 

ai = Ui x Uj_i , (33) 

Analytical work with the local winding axis is difficult, because 

a» • a i+ i = . l C ^\ 1 .' 1+ , 1 n ^ - cot(0 i+M ) cot(0i_i,i), (34) 

sin(6>i_i,i) sin(0 iii+ i) 

where we have defined cos(0j jJ ) = • Uj, giving rise to a four-point correlation function. The details 
of this result are shown in Appendix D. While directly computing the average of Eq. 34 is analytically 
intractable, the symmetry of the problem shows that, for a free WLC, (a, -aj+i) = (since Uj+i may be 
freely rotated about the Uj axis without changing 9 i:i+ i). The simulations show that the local winding 
axes for interior bonds are highly correlated for strongly confined chains, with (a 4 • a i+ i) collapsing on 
a single, increasing curve as a function of al p /R 2 (Fig. 7a). Correlations in the winding axis will thus 
develop more readily for smaller radii than will the oscillations seen in the bending correlation function 
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(depending on the ratio l p /R, see Eq. 32 and Fig. 6). The cndpoints of the chain are not strongly 
correlated to the interior behavior (Fig. 7a, inset), with a precipitous drop to (a.j • aj+i) < 0.1 at the 
cndpoints. This sharp drop suggests that the endpoints of the chain behave more like an unconfincd 
chain than do the interior monomers, with (a.j • a,+i) w 0, consistent with our physical picture of the 
origin of the endpoint effects (see the discussion above) . 

Correlations between the winding axes for interior monomers as a function of their separation As 
appear exponentially distributed (Fig. 7b), with a best fit 

(a( S )-a( S '))«(a l -a, +1 ) e -l As l/ 2 '°. (35) 

The sharp drop in (a(L/2) • a(s') near the cndpoints shows that the wrapping near the ends of the 
chain is uncorrelated to the interior wrapping, consistent with the behavior seen in the inset of Fig. 
7a. Forrey and Muthukumar [30] use (aj • z) as an order parameter in the study of the wrapping of 
DNA within the 4>29 phage (a natural choice, as the DNA loaded into the capsid along the z-axis). 
They find weak correlations between the local winding axis and the z-axis, as is expected due to the 
lack of correlations between the interior and the endpoints (Fig. 7) . 
Pressure Estimates: 

The free energy and pressure of the volume confined WLC can be computed using our mean field 
roots (Eq. 31). The exact expressions are somewhat lengthy, due to the endpoint terms involving po 
and p c , but in the limit of small R (relevant for most physical systems) we find 

/W-^ + t (W-^-f^+l, (36) 

with P = —dF/dV. The R^ 1 terms in the free energy and pressure are not present in the surface 
confined case, and are due entirely to the nonuniformity in (r 2 (s)) as a function of s (see Fig. 4), 
reflected in the fact that this term is proportional to p c . The coefficient of the term in Eq. 36 is 
negative, due to the fact that portions of the chain near the endpoints will be found on average closer 
to the center of the sphere than the interior monomers (as seen in Fig 4), resulting in a decrease in 
the pressure. The excellent agreement between simulations and the theoretical predictions (Fig. 8) 
shows that Eq. 36, with p rts 0.9, po ~ 0.95, and p c w 0.25 for strong confinement, can be used in the 
calculation of the entropy of confinement for a WLC. 

In order to determine the pressure directly from the simulations, we compute 

PA = J2^wairr (37) 

i 

with A = AttR 2 the surface area of the sphere, and fi^ wa ii the force of the i th monomer on the wall. 
In Fig. 8, we show the simulated results along with the full mean field expression for the pressure (of 
which Eq. 36 is the limit of small R). We find the agreement is excellent for a large range of L, l p and 
R, particularly for small R where endpoint effects are less important. 

We can compare our results to the experimental pressures determined by Smith et. al [25], using 
the 4>29 virus. The viral capsid is not spherical, with an icosohedral shell of radius w21nm and 
height «54nm, but has a volume equivalent to a sphere of radius ss 26nm. The fully packed virus 
contains a strand of DNA of length 6.6/im, with persistence length 50nm. If we neglect the excluded 
volume, electrostatic, and solvent-induced interactions of the DNA (a rather severe approximation), 
and take p s=s 0.9 (the saturating value of p, see Fig. 5a), we find P ss lkPa (=10 -3 pN/nm 2 ), 
almost 4 orders of magnitude lower than the 6MPa measured in the experiments. It is clear that the 
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behavior of a strongly confined wormlike chain is critically dependent on the intra-chain interactions, 
in agreement with a number of other studies [32, 33, 34, 27, 30, 38]. While we have found that entropy 
of confinement produces a negligible contribution to the experimentally observed pressure, excluded 
volume interactions will further restrict the conformational space available to the chain. Intra-chain 
interactions will lead to an increase in the cntropic contribution to the free energy and pressure, as 
has been seen in simulations [35]. Our results establish firmly, as noted some time ago [38], that the 
origin of spool- like order and the extremely large pressure of DNA in a capsid, is due to inter-segment 
and counterion-mediated interactions. 

5 Conclusions 

We have shown that WLCs in restricted spaces can be accurately treated by applying the mean field 
theory [44, 45, 46]. For a surface confined chain we can determine many average properties of the WLC, 
by replacing the rigid constraints of inextensibility (Ar 2 = a 2 ) and confinement (r 2 = R 2 ) with average 
constraints. We have shown that the mean field approach reproduces the exact results of Spakowitz 
and Wang [8], and reproduces the correct scaling coefficient of the free energy of confinement. The 
mean field approach is also able to determine the scaling and free energy of a surface confined WLC 
under tension, which may be of use in better understanding the wrapping of DNA around histones 
[2, 3]. The force-extension curve (FEC) for a strongly confined WLC differs greatly from the unconfincd 
FEC, with oscillatory behavior. 

We also find that the mean field method can approximately determine the behavior of a WLC 
confined to the interior of a sphere. Interior monomers (far from the endpoints) are approximately 
surface confined, with (r 2 (s)) w 0.9i? 2 for strongly confined chains, but endpoint effects dominate the 
behavior of the chain for s < 2R or s > L — 2R. Structural information about the confined chain can be 
determined by examining the correlations in the local winding axis, and we find that strongly confined 
stiff chains arc highly structured, even without intra-chain interactions. The mean field estimates of the 
pressure due to confinement show that the extreme pressures inside of a viral capsid are not strongly 
dependent on simple confinement entropy, but must arise from intra-chain and counterion-mediated 
interactions. The good agreement with the simulated pressures allows us to accurately estimate the 
free energy of confinement of a strongly confined WLC arising from energetic considerations alone as 

/JF«0.56^-l.l| + 31og(i2) (38) 

The excellent agreement between theory and simulations show that the mean field theory can be 
adapted to include thee effects of inter-segment interactions, even when semiflcxible chains are confined 
to restricted spaces. 

A The Q Matrix 

Defining x N = (x\, . . . ,Xn+i), we can rewrite the Hamiltonian in Eq. 2 as 

n ^ ' 
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where the elements of the tridiagonal matrix Q arc: 



Qi,i+2 — Qi+2,i 
Qt,»+1 — Qi+1,» 



lp 

2a^ 



2Xj 
a 2 
2Ai 



= 1,N 
else 



2a' 



+ a 2 ~r R 2 
A w , afcw+i 
a 2 + ^2 

2a 4 "T a 2 "T" fi 2 

a 4 "T a 2 _R 2 



2a 
5/„ 



i = 1 
i = iV+ 1 
i = 2,7V 
else 



(40) 
(41) 

(42) 



The structure of Q is unchanged under the transformation k\ = fcjv+i, = fcw, Ai = A^, ki = k for 
2 < i < N and X, = A for 1 < i < N. 

B Evaluation of the 1-D Confined Propagator 

We arc interested in evaluating the path integral in Eq. 6, 

x 2 (s) - X 



Z (x) = J V[x(a)] exp ^ - | J do x 2 {s) - X J ^ ds ±' \s) - J J ds:r 2 ( S )^, 



(43) 



subject to the boundary conditions x(0) = x , u(0) = u , x(L) = xl, and u{L) = ul- We write 
x(s) = f(s) + g(s), where g(0) = g(L) = g(0) = g(L) = and where 



f/ (4) (s)-A/( S ) + fc/( S ) = 0, 

with /(0) = Xo, f(L) = x L , /(0) = Mo, and f(L) = u L . If / satisfies the above differential equation, 
a simple integration by parts gives 



Z () (x) = /v(L).xp( -| 



- A[u L a; L - u a;o] 



tf(L) = Jv[g]e^^- l ^j\ S tf{ S )-xj\ S g 2 { S )-^j\ S g< 



(s) 



(44) 



(45) 



where g and g vanish at the boundaries. The exponential term in Eq. 44 can be evaluated by solving 
the differential equation for / directly, giving 



Z (x) = K(L)e 



•Mx 



where 



M= lpR 



2d(L) 



f Rmu Rm\2 Rvti\ 3 mu \ 

Rm\i Rrrin -m u — Rm\ 3 

Rmi3 —tnii ~m 33 /R m 3 ^jR 

\ -Rmi 3 m 3i /R m 33 /R J 



( -1 \ 

1 

-10 

y o i o o y 



(46) 



(47) 
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In Eq. 47, we have defined 

2u>iw 2 (l — cosh(LtJi) cosh(La;2)) + (ujf + w 2 ) sinh(Lwi) sinh(Lw2) 



d(L) = 

mu = 

m 12 = 

m 13 = 

mu = 

m 33 = 

m 34 = 



ui\ sinh(Lo>i) cosh(La;2) — ^2 cosh(Lwi) sinh(Lw 2 ) 

u> 2 sinh(Lw 2 ) ~ w i sinh(Lwi) 

, ,2 , ,2 

■ sinh(Lwi) sinh(Lw2) 



cosh(Lwi) — cosh(Loj 2 ) 
1 



1 



cosh(Lwi) sinh(Lcj 2 ) sinh(Lwi) cosh(.Lw2) 

— sinhfLwi) sinh(Lcj2), 

wi w 2 



with 



a;. 



1±A 1 



A 2 " 



(48) 



Note that the full propagator Z(x) in Eq. 10 can be written as Z{x) = Zq{x) exp(— x • Gx), with 
the matrix G containing terms suppressing excess endpoint fluctuations, 



/ li/R 2 l2 /R \ 

71 /R 2 -72 /R 

72/i? 8 

V -72/fl <y y 



(49) 



In general, computing average values involves calculating the determinant of M + G. Simplification 
of the determinant is a tedious process, but it is useful to note that 



Det(M) = A 2 - A 2 . 

.4, . mj, + + « -mum- 
^4 2 = 2mi 3 mi4 + mi 2 m 33 - mnm 3 4 



(50) 



with a similar relation holding for Dct(M + G). 

We can calculate K(L) by the evaluation of a simple integral. Following the standard method of 
Fcynman [49], we can write the propagator from (xo, uq) to (xl,Ul) as an integral over all intermediate 
points, (x s ,u s ), 



/oo 
dx s du s Z (x ,u ,x s ,u s ; s) x Z (x s ,u s ,x L ,u L ; L - s) 
-00 

= / / dx s du s K(s) exp I — • M(s)xi I 

J — 00 J — 00 V / 

x if (L - s) exp I - Xj • M(L - s) x 2 



(51) 
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where xi = (xo, x s , u , u s ) and x 2 = (x s ,Xl,u s ,Ul)- M has already been determined (Eq. 47), and 
K(L), given in Eq. 45, is independent of all xis and Ui's. The integral in Eq. 51 is tedious to evaluate, 
but yields 



Z (x;L) 



K(L)t 



•M(L) : 



K(s)K(L-s) 



2ir ( d(s)d(L-s) 
V\,(w?-w|)d(L) 



-x T -M(L)x 



We then find 



d(L) 



2 e- L ", 



(52) 



(53) 



where 77 is an arbitrary constant. 

In the limit of strong confinement (Ixjj 3> 1, see the main text) 



M 



lpRui\U)2 

2~ 



/ R(U)1 + LU 2 ) 



V 





R((Ji + u) 2 ) 




1 



(wf 1 + u)^)/R 



\ 

-1 





K cx V^ 2 (w 1 + w 2 )e- i ^ 1+W2 )/ 2 



(54) 



This strongly confined representation is significantly easier to work with when computing the mean 
field solutions. 

To ensure that our calculation of Z has the correct limiting behavior, we find 



lim 



dx L Z (x.) = e 



-kLz 



n 



2ir sin 



V 

h(ClL) J 



cxp 



2sinh(fiL) 



(55) 



(«o + u i) cosh(fiL) — 2uqUl 



+ 0(1/R)\, 



with fl = yf2\/l p , identical to the unconfincd propagator found in the work of Ha and Thirumalai 
[44, 46], except for the the term e - kLx l. Since k cx R~ 2 in all cases considered, the integration over 
£0 leads to a divergent integral. However, as R — > 00, the system becomes translationally invariant, 
so that integration over the initial position will be proportional to the radius of the confinement. The 
integral over the initial position then simply adds an irrelevant constant to the free energy, and we can 
write the one dimensional propagator 



1 f°° 

Zu(u ,u L )= lim — / dx dx L Z(x) 



(56) 



identical to the result found by Ha and Thirumalai [44, 46] up to a multiplicative constant. 

C Details of the simulations for volume confinement 

We have therefore performed a number of Langcvin Dynamics simulations with varying L, l p , and R 
when considering volume confinement. The Hamiltonian used is 



(5H 



k 
2^ 



JV j iV-l N+1 , 

^(|r i+ i - r;| - a) 2 - ^ ^ u 4 • u i+ i + es ^ ( 

i=0 i=0 i=0 ^ 



(R+a) 
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(57) 
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with L = Na. The first term ensures the connectivity of the chain, and we take k = 10 4 throughout, 
ensuring very stiff bonds. The second term accounts for the bending stiffness of the chain, with 
persistence length l p . We have confirmed directly that this Hamiltonian in the unconfined case (i.e. 
only the first two terms of Eq. 57 are used) gives (u(s) • u(s')) = e~l As l/ /p to within ~ 5%. The 
third term of the Hamiltonian approximately confines the chain to the interior of a sphere of radius 
R, using a Lennard- Jones repulsion. The confinement energy is on the order of es ksT when |rj| = R, 
and increases sharply for larger |rj|. We choose es = 1 throughout the simulations, which restricts 
|rj|/i? < 1.01 for all of the parameters we considered. To determine the equilibrium properties of the 
system, we use the low friction limit [52], with r\ — 0.1, and a timestep of h = 0.001 (in dimcnsionless 
units, or equivalently with the mass m = 1, spacing a = 1, and fc^T =1). In the simulations, we 
consider a chain with N = 200 for l p /a = 20, 50, and 100, withi?/a = 5, 6, 7, and 8. We also consider 
a chain with N = 100 for l p /a = 5, 10, 20, 50, and 100, with R/a = 5, 6, 7, 8, 9, 10, 11, 12, 15, 20, 
and 25. 



D Calculation of the winding axis 

Because of the spherical symmetry of the problem, we are free to choose our coordinate system such 
that it simplifies the calculation. We take Ui_i = z, defining the z-axis, and = (sin[Oi],O,cos[0i]), 
defining the cc-axis. We take our third bond to be Uj+i = (sin[6j + i] cos[y], sin[0 i+ i] sin[y>], cos[0,+i]). 
With 9ij the angle between bonds i and j, we see 0,-1,, = 6, and 8i-i.i+i = 0j+i. It is convenient to 
eliminate the azimuthal angle ip when computing cos(#^ + i) = • Uj+i, giving 

cos(<p) = cot(0 M+ i) cot(0j) - cos(0 i+ i) csc(6» M+ i) csc(0j) (58) 

In this coordinate system, = y, and |a i+ i| = sin(0j i j + i). To compute the dot product between the 
two winding axes, we need only a»+i • y = cos(0i) sin(0 i+ i) cos(<,o) — cos(0i + i) sin(0j). Eq. 34 is 
recovered upon substitution of cos(</?) in sn ■ a i+ i = a»+i • y/|a i+ i|. 

Acknowledgments: This work was supported in part by a grant from the National Science 
Foundation through NSF CHE 05-14056. 



17 



Figure Captions 

Fig 1: Representative structures for a WLC confined to the surface of a sphere of radius R = 3a. 
(a) shows l p = 2.5a and (b) shows l p = 20a. An enlargement of the polymer in (b) diagrams the 
positions and bond vectors Ui. 

Fig. 2: Free energy (3AF — f3F(R) — (3F(oo) as a function of R for a surface-confined WLC. /3F(oo) 
is determined from a simulation with R = 2 x 10 4 a. The symbols are the simulation data, where the 
lines are the theoretical results of eq. 15. Shown are Z p /a=20 (solid purple), 10 (dotted blue), 5 (dashed 
green), and 2.5 (dot-dashed red). The inset shows log(/3AF) as a function of R, displaying the good 
agreement between simulation and theory, particularly for large R. 

Fig. 3: Linear extension under an external tension of a surface confined WLC as a function of the 
radius. In (a), L = 450a and l p — 150a. The applied tensions are a(3f= 0.1 (solid blue), 0.05 (dashed 
green) and 0.01 (dotted red), displaying the oscillations in (xl — %o) for stiff chains. In (b), the same 
values of the force are applied to a chain of length L = 450a and l p = 15a, showing that the FEC of a 
flexible chain is monotonic. 

Fig. 4: (r 2 (s)) vs. s. (a): L = 100a and l p — 100a, with (from highest to lowest) R/a=5, 10, and 
20. The average monomer position is dominated by endpoint effects for s < 2R and s > L — 2R. (b): 
L = 200a and R — 5a, with (from highest to lowest) Z p /a=100, 50, and 20. With R fixed, variations 
in l p change only the value of p, but do not alter the behavior of the monomers. 

Fig. 5: (a) The mean field parameters p (solid blue), p (dotted red), and p c (dashed green) as a 
function of l p /R. Lines arc determined from a simulations with L = 100a for various l p and R. Points 
are from simulations with L — 200a. (b) a = l p /lo as a function of l p /R, determined by fitting eq. 32 
to the simulation results. Symbols are the fits for L = 100a, the line is the fit for L = 200a. 

Fig. 6: Bending correlation function (u(i/2) • u(s)) as a function of s for a chain with L = 200a. 
The points are simulation data, the solid lines are the theoretical results in Eq. 32. (a) has l p = 20a and 
R = 8a, with a w 2, and (b) has l p = 100a and R — 5 with a ~ 5/2. The agreement between theory 
and simulations is excellent, except near the endpoints. Representative configurations are shown to 
the right, and clearly displays the wrapping to the chain for a strongly confined WLC in (b). 

Fig. 7: (a) Average correlations in the nearest neighbor winding axis as a function of al p /R 2 , with 
the average taken over interior points only (i.e. 2R < s < L — 2R). The symbols are simulation 
data for L — 100a for various l p and R, with the line the simulation data for L = 200a. The inset 
shows the average nearest neighbor correlation for L = 200a, R = 5a, and l p /a = 100 (blue), 50 
(green), and 20 (red). The correlations drop sharply near the endpoints. (b) (a(L/2) • a(s)) as a 
function of s for L — 200a and R = 5a, for Z p /a=100, 50, and 20. The curves are a fit to the 
exponential (a(L/2) • a(s)) cx e~ L / 2l °. One representative configuration for L = 200a and l p = 100a, 
seen from two different viewpoints, clearly shows high correlations in the winding axis, except near the 
ends of the chain. 

Fig. 8: fiPV as a function of R for varying l p and L. In both, the dots are simulation results, and 
the solid line is the theoretical result, with p, po and p c taken directly from the simulation results, (a): 
L = 100a and l p /a~100 (red), 50 (green), 20 (blue), 10 (purple), and 5 (pink). The inset is a log-log 
plot for this data, (b): L = 200a and l p /a— 100 (red), 50 (green), and 20 (blue). 
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